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LATERAL BUCKLING OF CANTILEVERED I-BEAMS 
UNDER UNIFORM LOADS 


Stanley Poley, J. M. ASCE* 


INTRODUCTION 


In recent years, the engineer has become progressively more con- 
scious of the weight to rigidity ratio of structures. This has led him to 
reexamine questions of design of structural units from a more exact 
standpoint. As a result, he has encountered boundary value problems in 
elastic stability which are intractable by analytical means. Before the 
advent of the modern high speed computing machine, the labor involved 
in obtaining numerical solutions made numerical methods impractical 
in many important cases. An example of one of these is the problem of 
lateral buckling of a cantilevered I-beam subjected to a uniform load. 
A numerical solution of this problem, obtained with the aid of the IBM 
Card Programmed Calculator, is presented herein. 


Differential Equations and Boundary Conditions 


With coordinate axes as shown in Figure 1, the total energy of an 
I-beam buckled laterally under the action of a uniform normal load act- 
ing along the centroidal axis of the cross section is given by (1) 


where 
= minimum flexural rigidity of beam. 
U = lateral deflection of beam in plane of minimum rigidity. 
L = length of beam. 
D = flexural rigidity of one flange in its own plane. 
h = depth of beam. 
f = angular rotation of section of abscissa Z . 
C = torsional rigidity of cross-section. 
M = ¥ (i-Z)* = bending moment in the plane of maximum flexural 
rigidity of section of abscissa Z . 
= intensity of uniform normal load. 

Eq. (1) is derived under the same assumptions used by Timoshenko 
(2) in his “Theory of Elastic Stability”, Chapter V; it applies therefore 
to doubly symmetrical I-sections with the principal moment of inertia 
I y much smaller than I x. 
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Applying the usual methods of the calculus of variations to Eq. (1), 
we may derive the differential equation governing B ° 


w mM? =O 
and the boundary conditions 


in Eq. (2), the buckling equation becomes:* | 


while the boundary conditions (3) are: 
(to) = = = ae 


Solution by Finite Differences 


Subdividing the interval [0,1] into Y) equal parts of length h= Yn and 
taking the first terms in the expansions of the derivatives in terms of 
central differences (3), Eq. (5) becomes: 


- h*ti-x:)* pi=0 


where (i = A(%), Xi being a pivotal abscissa of the interval [0,1]. Substi- 
tuting in Eq. (7) the values of the central differences in terms of pivotal 
values, we obtain: 


SPs + [2€5-1)- — + Pisa =O (8) 


3. From here on, primes indicate differentiation with respect to X. 
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Similarly, the boundary conditions (8) become: 
| =P 
= ZPn 
= 2(S-2)( Bn-fn-1) + 


From equations (4a), (4b) and (9b), the critical total load takes the 
form 


(11) 


(12) 


To solve the boundary value problem defined by Eq’s. (8), (10) one 
fixes a value of */g?, and a value of N. Eq. (8) is then applied at the 
pivotal X; ({=1,2,----,M) using Eq. (10) to eliminate B,, 2, Aner, Bnez - 
This yields a set of Nl linear, algebraic, homogeneous equatigns in the 
Vi unknowns f,-----,@n and this system has solutions A, not identically 
zero, if and only if its determinant A , which is a function of K’, is 
equal to zero. In problems of elastic stability, the designer will be 
mainly interested in the lowest buckling load, and hence the determinen- 
tal equation 


A(k’) =O (13) 


is solved for its smallest root. Having done this for one value of N, the 
operation is repeated with a larger Y) to obtain a better approximation to 
the true value of K*. By increase in the number N of pivotal points, K? 
may be evaluated, theoretically, to any degree of accuracy. 

For large values of N , the evaluation of the roots of A(K*) =O is im- 
practical without the use of special computing devices. The actual cal- 
culations for this problem were performed on the LB.M. Card Pro- 
grammed Electronic Calculator, Model II, which was made available to 
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where 
$=4+(th) (9) 
out 
K*= (h) (9b) 
(10) 
| 
where 
m= 
(t/a) 


the author by the Watson Scientific Computing Laboratory, at Columbia 
University. A detailed account of the method employed for these calcu- 
lations is to be found else where (4), 


Results 


The following table gives the factor M versust’/a*. For a given I- 
beam, 1?/a? is computed by means of the formula 


The corresponding value of ™ is obtained from the table, and the criti- 
cal total load is calculated from Eq. (11). 


m [66.9 [50.5 [440] 99.0 [54.9] 31.9 [20.6 
iz [16 [32 | 40 [100 | 
m_[2e.3 [27.1 [a6.o [29.5 [21.9 [20.9 [17.6 | 


We see that, as the ratiot*/o increases, the factor ™ diminishes and ap- 
proaches as a limit the value ™ = 12.85, obtained for lateral buckling of 
a cantilivered rectangular beam under similar loading conditions (5), 


Example 


Consider a 24I79.9 section, with L = 20 ft. The flexural rigidity D 
of a flange can be taken equal to half of the flexural rigidity of the I- 
beam in the lateral direction, while the torsional rigidity ( may be cal- 
culated from the approximate formula (6), 


(abt? +ht, ) 


where b is the flange width, t the average flange thickness, h the beam 
depth, t, the web thickness, and G the shear modulus. For the above 
section, b= 7", t = 0.87", t,= 0.5", h= 24", B = 42.9E in4# and D= 21.5E 
in4, E being the elastic modulus. Taking E = 2.66, we get C # 1.57 in4. 
Thus 


DF 21.5+2* 


From the table, we obtain M* 26.8. With E = 30 x 106 psi, we have 


(gh = = 20.8 ™ 


(20«12)? 


a 
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Figure 1 - Coordinate Axes 
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